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1.  INTRODUCTION 


The  U.S.  Army  Research  Laboratory  (ARL)  completed  an  initial  investigation  of  the  flow  field  within 
a  typical  test  chamber  operated  by  the  Army  Edgewood  Research,  Development  and  Engineering  Center 
(ERDEC).  The  ERDEC  test  chamber  is  designed  to  mix  compressor-driven  airflow  with  gas/solid  effluent 
from  a  test  article  placed  inside  the  chamber.  An  example  of  such  a  test  article  is  a  smoke  generator,  or 
smoke  pot,  commonly  used  on  the  battlefield  to  provide  a  means  of  obscurant  During  the  test,  the 
air/effluent  flow  field  is  exhausted  from  the  test  chamber  for  analysis.  In  order  to  simulate  this  flow,  the 
ARL  applied  computational  fluid  dynamics  (CFD)  codes  that  include  multispecies  chemical  kinetics  as 
well  as  multiphase  (particulate)  submodels.  These  codes  were  developed  at  ARL  to  numerically  solve  the 
Navier-Stokes  equations  and  simulate  the  chemically  reacting,  multiphase  flow  field  in  gun  propulsion 
systems.  This  code  has  been  used  successfully  for  other  applications  at  ARL  (Nusca  1989,  1991, 1993). 

Application  of  the  code  to  the  present  study  involved  generating  a  computational  mesh  that  covered 
the  chamber  interior  as  well  as  specifying  proper  boundary  conditions  on  the  chamber  walls,  chamber  top 
(air  inflow),  chamber  exit  (outflow),  and  test  article  (effluent  outflow),  as  depicted  in  Figure  1.  The 
governing  equations,  boundary  conditions,  and  solution  method  are  outlined  in  this  report.  Numerical 
solutions  of  the  gas  flow  and  effluent  concentration  distributions  in  the  test  chamber  were  generated  for 
operating  times  up  to  4.5  min.  Graphical  results  with  discussion  are  presented  in  this  report.  Numerical 
simulations  reveal  that  certain  values  of  chamber  through-flow  induce  flow  patterns  within  the  chamber 
that  are  dominated  by  vortices.  This  flow  pattern  increases  the  effluent  residence  time  in  the  chamber  as 
well  as  the  mixing  of  gas/particulate  from  the  test  article  with  air.  The  test  article  effluent  jet  feeds 
effluent  into  this  vortical  motion,  and  only  that  flow  that  is  trapped  near  the  chamber  floor  is  drawn  out 
of  the  chamber.  Pockets  of  high  effluent  concentration  can  form  in  the  chamber. 

2.  GOVERNING  EQUATIONS 

For  purposes  of  producing  a  timely  initial  investigation,  the  cylindrical  test  chamber  was  modeled  as 
two-dimensional  (2D).  The  governing  equations  are  written  in  Cartesian  coordinates  with  velocity 
components  u  and  v  for  the  x  (along  chamber  floor)  and  y  (along  chamber  height)  directions,  respectively 
(see  Figure  1).  The  Reynolds-Averaged  Navier-Stokes  (RANS)  equations  describe  the  2D  reacting  gas 
flow  (N  species  mixture)  in  the  chamber  given  conditions  at  the  boundaries  of  the  geometry.  These  partial 
differential  equations  describe  the  time  (t)  evolution  of  the  dependent  variables  of  velocity  (u,  v),  pressure 
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(p),  mixture  density  (p),  species  mass  fraction  (Oj,  for  i  =  1  to  N  species),  internal  energy  (e),  temperature 
(T,  derived  from  energy),  and  viscous  shear  stresses  (t). 
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The  shear  stress  terms  are  given  by 
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In  these  equations,  Oj  and  (Oj  are  the  mass  fraction  and  chemical  production  terms  for  the  i*  species.  For 
the  present  application,  finite-rate  chemical  production  terms  were  not  used.  Chemical  reaction  was 
modeled  as  an  infinitely  fast,  one-step,  unidirectional  (i.e.,  forward)  reaction  of  smoke  pot  effluent  (i  =  1) 
and  air  (i  =  2)  to  form  product  (i  =  3)  for  stoichiometric  air/effluent  ratio  of  0.17  and  effluent  density 
above  50  g/m^.  The  reaction  temperature  was  taken  as  680”  C. 
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The  temperature  dependence  of  the  species  viscosity,  pj,  and  thermal  conductivity,  iq,  can  be  modeled 
using  Sutherland’s  law  (White  1974), 
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The  terms  p^,  T^,  and  S  can  vary  with  species  but  were  assumed  to  be  constant  with  values  of  Sj,  =  199  R, 
Top  =  491.6  R,  Po  =  0.1716  mP,  =  350  R,  T^^  =  491.6  R,  =  0.0139  BTU/h-ft-R.  The  mixture 
viscosity  and  thermal  conductivity  (mixture  quantities  are  denoted  by  subscript  m)  are  determined  using 
WUke’s  law  (Wilke  1950),  denoting  f  as  p  or  k. 
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where  Xj  and  M^  are  the  mole  fraction  (Xj  =  po/Mj)  and  molecular  weight  of  the  i*^  species,  respectively 
(Ml  =  97.94,  M2  =  28.8,  and  M3  =  63.37  g/mole).  Pick’s  law  (White  1974)  is  used  to  relate  the  mixture 
diffusivity  to  the  mixture  viscosity  through  the  Schmidt  number,  Sc  =  Pjj/(p  D),  assumed  unity.  The 
specific  heat  at  constant  pressure  of  each  species  (per  mass)  is  generally  given  by  the  following  fourth- 
order  polynomial  curve  fit  (Drummond,  Rogers,  and  Hussaini  1987): 

Q 

_ZL  =  A:  +  BjT  +  CT^  +  DjT^  +  EjT^. 

Rill  1 


For  the  present  study,  Cp  was  assumed  constant  with  values  Cpi  =  0.2878,  Cp2  =  0.238,  and 
Cp3  =  0.1277  cal/g°C.  The  mixture  pressure  (equation  of  state),  enthalpy,  total  energy  per  unit  volume. 
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and  ratio  of  specific  heats  are  given  by  (R^  is  the  universal  gas  constant  and  AHg  is  the  heat  of  formation 
for  species  i) 
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An  algebraic  turbulence  model  (Bradshaw,  Cebeci,  and  Whitelaw  1981)  was  used.  In  this  model,  the 
eddy  viscosity,  Pj,  is  computed  assuming  that  the  viscous  layer  consists  of  an  inner  and  an  outer 
component.  The  inner  region  follows  the  Prandtl  mixing  length  formulation  based  on  a  prescribed 
characteristic  length  scale,  L,  a  boundary  layer  intermittency  factor,  e  (having  a  value  of  0  for  laminar, 
1  for  turbulent  flows,  and  a  function  of  x  for  transitional  flows),  the  displacement  thickness  of  the  layer, 
5,  and  a  constant,  a. 

(pt).  =  L^y  ||i^i  II ,  Oiyiiyc  . 
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Here,  y^,  is  a  prescribed,  small  distance  from  the  solid  boundary,  and  y^  is  the  edge  of  the  viscous  layer. 
Further  details  can  be  obtained  from  Bradshaw,  Cebeci,  and  Whitelaw  (1981).  The  fluid  viscosity  is  then 
p  =  Pf  where  obtained  using  Sutherland’s  law  and  Wilke’s  law. 

3.  BOUNDARY  CONDITIONS  AND  INITIAL  CONDITIONS 

The  boundaries  of  the  test  chamber  (see  Figure  1)  are  the  air  inlet  at  the  top  (roof),  the  exit  port  on 
the  chamber  floor  (connected  by  ducts  to  the  wind  tunnel  fan),  and  the  vertical  walls.  The  smoke  pot  is 
placed  on  the  chamber  floor,  near  the  chamber  exit  port.  Since  the  governing  equations  are  elliptic  (low- 
speed  flow),  conditions  along  these  boundaries  must  prescribe  values  of  the  dependent  variables,  the 
gradient  of  the  dependent  variables  in  the  boundary-nonnal  direction,  or  an  algebraic  relation  which 
connects  the  values  of  the  dependent  variables  to  the  normal  component  of  velocity. 

At  the  air  inlet,  x-direction  profiles  of  all  dependent  variables,  p,  u,  v,  a,  T,  and  p,  are  specified.  It 
is  assumed  that  the  flow  at  the  inlet  consists  of  air  and  that  convection/diffusion  of  effluent  to  the  chamber 
top  is  not  permitted  to  exit  the  chamber.  By  mass  conservation,  the  inlet  flow  velocity  was  specified  as 
u  =  .062  ft/s,  and  a  parabolic-shaped  profile  was  assumed. 

The  exit  port  velocity  was  specified  as  u  =  2.96  ft/s  (5,380  1/min)  with  a  parabolic-shaped  profile. 
Boundary-normal  gradients  of  all  dependent  variables  at  the  exit  plane  are  zero.  Mass  that  exits  the  port 
is  not  assumed  to  reenter. 

The  no-slip/no-penetration  condition  (u  =  v  =  0)  is  applied  to  the  solid  chamber  and  smoke  pot  walls. 
The  walls  are  assumed  to  be  adiabatic  (i.e,  normal  derivative  of  T  set  to  zero).  The  normal  gradient  of 
all  mass  fractions,  da^dn,  are  also  set  to  zero. 

The  top  of  the  smoke  pot  was  assumed  to  be  a  constant  mass  flux  source  of  effluent  with  u  =  12.7  ft/s, 
T  =  320°  C,  Ml  =  97.94  g/mole,  and  Cp  =  0.2878  cal/g°  C. 

4.  COMPUTATIONAL  ALGORITHM 

Equation  (1)  can  be  reduced  to  a  successive-substitution  formula  for  a  general  dependent  variable,  W, 
at  each  node  on  the  computational  grid.  Central  finite-differences  are  used  for  the  diffusive  (arrays  Gj 
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and  Gj)  and  source  terms  (array  Q)  and  upwind  differences  for  the  convective  terms  (arrays  Fj  and  Fj). 
Using  upwind  differencing  in  the  species  conservation  equations  (i.e.,  W  =  pOj)  reduees  the  occurrence 
of  negative  species  mass  fractions  in  mixing  layers.  The  resulting  system  of  equations  for  the  entire  grid 
is  solved  using  a  Gauss-Seidel  relaxation  scheme.  Each  iteration  cycle  is  made  up  of  J  subcycles,  where 
J  is  the  number  of  equations  being  considered.  In  each  subcycle,  grid  points  are  scanned  row  by  row,  and 
a  single  variable  is  updated.  When  aU  subcycles  are  completed,  a  new  iteration  cycle  in  which  the  values 
of  the  variables  from  the  latest  iteration  are  immediately  used  is  started.  This  is  consistent  with  the 
Gauss-Seidel  methodology.  Convergence  is  satisfied  when  the  greatest  relative  change  in  any  flow 
variable  is  smaller  than  a  prescribed  tolerance.  See  Nusca  (1989,  1991)  for  further  details. 

5.  RESULTS  AND  DISCUSSION 

Figure  1  shows  the  computational  grid  used  to  discretize  the  chamber  interior.  The  number  of  grid 
nodes  in  the  x  and  y  directions  are  75  and  50,  respectively  (3,750  nodes  total).  Grid  node  clustering  was 
used  to  resolve  flow  gradients  near  the  smoke  pot. 

The  simulation  was  run  for  approximately  1  min  to  establish  steady  flow  in  the  chamber  before  the 
smoke  pot  was  activated.  Figure  2  shows  the  streamline  (contour  lines  of  constant  stream  function) 
patterns.  Note  that  a  large  counterclockwise  vortex  resides  to  the  upper  left  of  the  smoke  pot  (established 
by  flow  from  the  chamber  inlet  that  must  turn  at  the  chamber  floor)  and  that  a  smaller  clockwise  vortex 
resides  over  the  smoke  pot  (established  by  flow  rising  in  the  vertical  direction  that  is  turned  by  the 
chamber  iidet  flow  at  the  top). 

Figures  3,  4,  5,  and  6  show  the  flow  streamline  pattern  after  1,  2,  3,  and  4  min  of  smoke  pot 
operation,  respectively.  Initially,  flow  from  the  smoke  pot  rises  toward  the  chamber  top,  establishing  two 
small  vortices  near  the  pot,  rotating  in  opposite  directions.  At  later  times,  the  flow  settles  into  a  large 
counterclockwise  vortex  offset  from  the  centerline  of  the  chamber  and  fed  by  the  smoke  pot  jet  Flow 
entrained  in  the  chamber  exit  port  is  limited  to  that  trapped  near  the  chamber  floor.  Figure  7  shows  the 
flow  streamline  pattern  at  4.5  min,  which  is  0.5  min  after  the  smoke  pot  has  ceased  operation.  The  vortex 
has  reduced  in  size  and  is  centered  between  the  vertical  chamber  walls. 

Figures  8-12  and  Figures  13-17  show  contours  of  smoke  pot  effluent  mass  fraction,  Oi,  (mass  of 
effluent/total  mass)  and  effluent  density  (product  of  mass  fraction  and  mixture  density),  respectively,  at 
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times  1-4.5  min.  At  early  times,  effluent  concentrations  are  high  in  the  smoke  pot  jet.  At  later  times, 
the  effluent  is  entrained  in  the  chamber  vortex  and  diffused  to  smaller  concentrations.  Even  at  later  times, 
pockets  of  high  concentration  (50  gAn^  or  greater)  can  be  noted.  The  flow  pattern  is  not  greatly  disturbed 
by  the  chamber  exit  port  on  the  floor.  Figures  12  and  17  show  the  effluent  mass  fraction  and  density  at 
4.5  min,  0.5  min  after  flow  from  the  smoke  pot  has  been  stopped.  The  chamber  vortex  has  swept  effluent 
into  the  vicinity  of  the  smoke  pot  where  it  becomes  trapped  at  large  concentration  levels. 

6.  CONCLUSIONS 

Due  to  the  low-vertical  flow  velocity  (0.06-2.7  ft/s)  through  the  chamber  induced  by  the  small 
chamber  exit  port  on  the  floor,  the  natural  flow  pattern  in  the  chamber  is  one  that  is  dominated  by  rotating 
vortices.  This  pattern  increases  the  flow  residence  time  in  the  chamber  and  mixes  gases  from  the  smoke 
pot  with  air  (similar  to  a  "weU-stirred  reactor").  The  smoke  pot  jet  feeds  effluent  into  this  vortical  motion, 
with  only  that  flow  that  is  trapped  near  the  chamber  floor  exiting  the  chamber.  As  a  result,  effluent  is 
allowed  to  fonn  pockets  of  high  concentration  that  may  chemically  react  with  the  fresh-air  supply  fed 
from  the  chamber  inlet  (i.e,  top).  After  the  smoke  pot  ceases  operation,  the  chamber  vortex  concentrates 
effluent  near  the  chamber  wall.  A  larger  chamber  exit  port  and  forced  exit  velocity  (controlled  by  the 
wind  tunnel  fan)  may  assist  in  breaking  these  vortices  and  evacuating  the  chamber  at  the  higher  rate.  The 
increased  chamber  through-flow  should  be  sufficient  to  turn  the  smoke  pot  jet  toward  the  exit.  Numerical 
simulations  aimed  at  predicting  this  effect  have  not  been  pursued. 

The  numerical  simulations,  results,  discussions,  and  conclusions  reached  in  this  report  are  subject  to 
the  assumptions  used  in  the  model  and  the  information  supplied  to  the  model  in  the  form  of  boundary 
conditions.  While  the  confidence  level  in  the  model  is  high  (based  on  performance  in  simulating  other 
problems),  further  studies  that  test  model  sensitivity  to  the  supplied  boundary  conditions  should  be 
conducted.  A  full  three-dimensional  simulation  should  be  conducted  to  model  the  perforated  chamber  top 
wall,  in  addition  to  three  discrete  smoke  pot  exit  holes  as  well  as  flow  obstructions  (i.e.,  pipes)  in  the 
chamber.  These  are  thought  to  represent  secondary  effects  in  the  simulation  of  unknown  final  effect  on 
the  results. 
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Figure  1.  Schematic  of  test  chamber  and  smoke  pot  showing  computational  grid. 
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Figure  2.  Stream  function  contours  before  smoke  pot  operation. 
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Figure  3.  Stream  function  contours  1,0  min  after  smoke  pot  startup. 
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Figure  6.  Stream  fiinction  contours  4.0  min  after  smoke  pot  startup. 
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Figure  7.  Stream  function  contours  4.5  min  after  smoke  pot  startup.  0.5  min  after  smoke  pot  termination. 
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Figure  8.  Effluent  mass  fraction  contours  1.0  min  after  smoke  pot  startui 
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Figure  9.  Effluent  mass  fraction  contours  2.0  min  after  smoke  pot  staitui 
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Figure  10.  Effluent  mass  fraction  contours  3.0  min  after  smoke  pot  startui 
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Figure  11.  Effluent  mass  fraction  contours  4.0  min  after  smoke  not  startun. 
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Figure  13.  Effluent  density  (g/m^)  contours  1.0  min  after  smoke  pot  startui 
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Figure  14.  Effluent  density  (gAn-^l  contours  2.0  min  after  smoke  pot  startui 
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Figure  15.  Effluent  density  (e/mh  contours  3.0  min  after  smoke  pot  startun. 
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Figure  17.  Effluent  density  contours  4.5  min  after  smoke  pot  startup.  0.5  min  after  smoke  pot  termination. 
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LIST  OF  SYMBOLS 


Cp  =  specific  heat  capacity,  constant  p 

Cy  =  specific  heat  capacity,  constant  volume 

D  =  mass  diffusion  coefficient 

e  =  specific  total  internal  energy 

F,  G  =  flux  vectors 

h  =  molar  specific  enthalpy 

L  =  Prandtl  mixing  length 

M  =  molecular  weight 

N  =  total  number  of  species 

p  =  static  pressure 

R  =  specific  gas  constant,  (y-l)c^y 

R„  =  universal  gas  constant,  R 

Sc  =  Schmidt  Number,  D 

t  =  time 

T  =  static  temperature 

u  =  axial  velocity 

V  =  radial  velocity 

W  =  dependent  variable  vector 

X,  y  =  Cartesian  coordinates 

X  =  species  mole  fraction 


Greek  S5 


y  =  ratio  of  specific  heats,  Cp/c^ 

AHf  =  enthalpy  of  formation 

5  =  boundary  layer  displacement  thickness 

e  =  boundary  layer  intermittency  factor 

K  =  heat  transfer  coefficient 

p  =  molecular  viscosity 
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p  =  density 

o  =  species  mass  fraction 

X  =  shear  stress  tensor 

0  =  chemical  production  term 

=  source  term  vector 


e  =  edge  of  the  viscous  layer 

i  =  i*  species 

m  =  mixture  quantity 

p  =  constant  pressure 

t  =  turbulence  quantity 

V  =  constant  volume 

X  =  x-direction 

y  =  y-direction 
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5 

DIRUSARL 

ATTN  AMSRL  OP  AP  L 

1 

DIR  USA  ADVANCED  SYSTEMS 

ATTN  AMSAT  R  NR  MS  219  1 

R&A  OFC 

AMES  RESEARCH  CENTER 

MOFFETT  FLD  CA  94035-1000 
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1  HQDA  (SARD-TR/MS.  K.  KOMINOS) 

WASH  DC  20310-0103 

1  HQDA  (SARD-TR/DR.  R.  CHATT) 

WASH  DC  20310-0103 

1  SDOP/TNI 

ATTN:  L.  H.  CAVNEY 
PENTAGON 

WASHINGTON  DC  20301-7100 

6  COMMANDER,  U.S.  ARMY  ARDEC 

ATTN:  SMCAR-AET-A, 

R.  DEKLEINE 
R.  KLINE 

R.  BOTTICELLI 
H.  HUDGINS 

J.  GRAU 

S.  KAHN 

PICATINNY  ARSENAL  NJ  07806-5001 

1  COMMANDER,  U.S.  ARMY  ARDEC 
ATTN;  SMCAR-CCH-V,  PAUL  VALENTI 
PICATINNY  ARSENAL  NJ  07806-5001 

2  COMMANDER 

U.S.  ARMY  RESEARCH  OFFICE 
ATTN:  TECHNICAL  LIBRARY 
D.  MANN 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC  27709-2211 

1  DIRECTOR 

U.S.  ARMY  RESEARCH  OFFICE 
ATTN:  AMXRO-MCS,  MR.  K.  CLARK 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC  27709-2211 

1  DIRECTOR 

U.S.  ARMY  RESEARCH  OFFICE 

ATTN:  AMXRO-RT-P,  LIBRARY  SERVICES 

P.O.  BOX  1211 

RESEARCH  TRIANGLE  PARK  NC  27709-2211 

4  COMMANDER 

NAVAL  RESEARCH  LABORATORY 
ATTN:  TECHNICAL  LIBRARY 
CODE  4410, 

K.  KAILASANATH 
J.  BORIS 

E.  ORAN 

WASHINGTON  DC  20375-5000 


1  COMMANDER 

NAVAL  SURFACE  WARFARE  CENTER 
ATTN:  DR.  F.  MOORE 
DAHLGRENVA  22448 

7  COMMANDER 

NAVAL  SURFACE  WARFARE  CENTER 
ATTN:  T.  C.  SMITH 
K.  RICE 
S.  MITCHELL 
S.  PETERS 
J.  CONSAGA 

C.  GOTZMER 
TECHNICAL  LIBRARY 

INDIAN  HEAD  MD  20640-5000 

2  COMMANDER 

NAVAL  SURFACE  WARFARE  CENTER 
ATTN:  CODE  R44,  DR.  A.  WARDLAW 
K24,  B402-12,  DR.  W.  YANTA 
WIDTE  OAK  LABORATORY 
SILVER  SPRING  MD  20903-5000 

1  USAF  WRIGHT  AERONAUTICAL 

LABORATORIES 

ATTN:  AFWAL/FMG,  DR.  J.  SHANG 
WPAFBOH  45433-6553 

1  AFOSR/NA 

ATTN:  J.  TISHKOFF 
BOLLING  AFB  DC  20332-6448 

3  AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/FXA, 

STEPHEN  C.  KORN 
BRUCE  SIMPSON 
DAVEBELK 

EGUNAFBFL  32542-5434 

2  WL/MNSH 

ATTN:  R.DRABCZUK 

D.  LITTRELL 

EGLINAFBFL  32542-5434 

1  LOS  ALAMOS  NATIONAL  LABORATORY 

ATTN:  MR.  BILL  HOGAN 
MS  G770 

LOS  ALAMOS  NM  87545 
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2  DIRECTOR 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  DIV.  1554,  DR.  W.  OBERKAMPF 
DIV.  1554,  DR.  F.  BLOTTNER 
ALBUQUERQUE  NM  87185 

4  DIRECTOR 

NASA,  LANGLEY  RESEARCH  CENTER 
ATTN:  TECH  LIBRARY 

MR.  D.  M.  BUSHNELL 
DR.  M.  J.  HEMSCH 
DR.  J.  SOUTH 
LANGLEY  STATION 
HAMPTON  VA  23665 

2  DIRECTOR 

NASA,  LANGLEY  RESEARCH  CENTER 
ATTN:  MS  408, 

W.  SCALLION 
D.  WITCOFSKI 
HAMPTON  VA  23605 

6  DIRECTOR 

NASA,  AMES  RESEARCH  CENTER 
ATTN:  MS-227-8,  L.  SCHIFF 
MS-258-1,  T.  HOLST 
MS-258-1,  D.  CHAUSSEE 
MS-258-1,  M.  RAI 
MS-258-1,  P.  KUTLER 
MS-258-1,  P.  BUNING 
MOFFETT  FIELD  CA  94035 

1  UNITED  STATES  MILITARY  ACADEMY 

DEPARTMENT  OF  MECHANICS 
ATTN:  LTC  ANDREW  L.  DULL 
WEST  POINT  NY  10996 

1  UNIVERSITY  OF  CALIFORNIA,  DAVIS 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN:  PROF.  H.  A.  DWYER 
DAVIS  CA  95616 

3  SCIENCE  AND  TECHNOLOGY  INC. 
ATTN:  DR.  ALAN  GLASSER 

MR.  BRUCE  LOHMAN 
MR.  DAVE  MAURm 
4001  NORTH  FAIRFAX  DR.,  NO.  700 
ARLINGTON  VA  22203-1618 


1  MASSACHUSETTS  INSTITUTE  OF 

TECHNOLOGY 
ATTN:  TECH  LIBRARY 
77  MASSACHUSETTS  AVE. 

CAMBRIDGE  MA  02139 

1  GRUMANN  AEROSPACE  CORPORATION 

AEROPHYSICS  RESEARCH  DEPARTMENT 
ATTN:  DR.  R.  E.  MELNIK 
BETHPAGENY  11714 

1  AEDC 

CALSPAN  FIELD  SERVICE 
ATTN:  MS  600,  DR.  JOHN  BENEK 
TULLAHOMATN  37389 

1  VIRGINIA  POLYTECHNIC  INSTITUTE 

AND  STATE  UNIVERSITY 
DEPARTMENT  OF  AEROSPACE  AND  OCEAN 
ENGINEERBSIG 

ATTN:  DR.  CLARK  H.  LEWIS 
BLACKSBURG  VA  24061 

1  ADVANCED  TECHNOLOGY  CENTER 

ARVIN/CALSPAN 

AERODYNAMICS  RESEARCH  DEPARTMENT 
ATTN:  DR.  M.  S.  HOLDEN 
P.O.  BOX  400 
BUFFALO  NY  14225 

1  THE  PENNSYLVANIA  STATE  UNIVERSITY 

DEPT  OF  AEROSPACE  ENGINEERING 
ATTN:  DR.  G.  S.  DULKRAVICH 
UNIVERSITY  PARK  PA  16802 

3  THE  PENNSYLVANIA  STATE  UNIVERSITY 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN:  V.  YANG 
K.  KUO 
C.MERKLE 

UNIVERSITY  PARK  PA  16802-7501 

1  UNIVERSITY  OF  ILLINOIS  AT  URBANA 

CHAMPAIGN 

DEPARTMENT  OF  MECHANICAL  AND 
INDUSTRIAL  ENGINEERING 
ATTN:  DR.  J.  C.  DUTTON 
URBANA  IL  61801 
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1  UNIVERSITY  OF  MARYLAND 

DEPT  OF  AEROSPACE  ENGINEERING 
ATTN:  DR.  J.  D.  ANDERSON,  JR. 
COLLEGE  PARK  MD  20742 

1  UNIVERSITY  OF  NOTRE  DAME 

DEPT  OF  AERONAUTICAL  AND 
MECHANICAL  ENGINEERING 
ATTN:  PROF.  T.  J.  MUELLER 
NOTRE  DAME  IN  46556 

1  UNIVERSITY  OF  TEXAS 

DEPT  OF  AEROSPACE  ENGINEERING 
MECHANICS 

ATTN:  DR.  D.  S.  DOLLING 
AUSTIN  TX  78712-1055 

1  UNIVERSITY  OF  DELAWARE 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN:  DR.  JOHN  MEAKIN,  CHAIRMAN 
NEWARK  DE  19716 

1  ARROW  TECHNOLOGY  ASSOCIATES,  INC. 

ATTN:  W.  HATHAWAY 
P.O.  BOX  4218 

SOUTH  BURLINGTON  VT  05401-0042 

1  PAUL  GOUGH  ASSOCIATES,  INC. 

ATTN:  P.S.  GOUGH 

1048  SOUTH  ST. 

PORTSMOUTH  NH  03801-5423 

2  PRINCETON  COMBUSTION  RESEARCH 

LABORATORIES,  INC. 

ATTN:  N.  MER 

N.  A.  MESSINA 

PRINCETON  CORPORATE  PLAZA 
11  DEERPARK  DR.,  BLDG.  IV,  SUITE  119 
MONMOUTH  JUNCTION  NJ  08852 

1  ROCKWELL  INTERNATIONAL 

ROCKETDYNE  DIVISION 
ATTN:  BA08,  R.  B.  EDELMAN 
6633  CANOGA  AVE. 

CANOGAPARKCA  91303-2703 


1  VERTTAY  TECHNOLOGY,  INC. 
ATTN:  E.  nSHER 
4845  MILLERSPORT  HIGHWAY 
EAST  AMHEREST  NY  14501-0305 


ABERDEEN  PROVING  GROUND 
31  DIR,  USARL 

ATTN:  AMSRL-WT-P,  MR.  ALBERT  HORST 
AMSRL-WT-PB, 

DR.  E.  SCHMIDT 
MR.  B.  GUIDOS 
DR.  P.  PLOSTINS 
DR.  J.  SAHU 
MR.  P.  WEINACHT 
DR.  G.  COOPER 
AMSRL-WT,  DR.  A.  BARROWS 
AMSRL-WT-PD,  DR.  B.  BURNS 
AMSRL-WT-PA, 

DR.  T.  MINOR 

MR.  M.  NUSCA  (7  CP) 

MS.  G.  WREN 
DR.  T.  COFFEE 
MR.  J.  DESPIRITO 
DR.  D.  KOOKER 

MR.  D.  KRUCZYNSKI 
DR.  G.  KELLER 

DR.  F.  LBERATORE 
AMSRL-WT-W,  DR.  C.  MURPHY 
AMSRL-WT-WB,  DR.  W.  D’AMICO 
AMSRL-WT-NC, 

MS.  D.  HISLEY 
MR.  R.  LOTTERO 
DR.  K.  OPALKA 

AMSRL-CI-C, 

DR.  W.  STUREK 
DR.  N.  PATEL 

3  DIRUSAERDEC 

ATTN  SCBRD-RTB, 

MR.  D.  PALUGHI 
MR.  L.  BICKFORD 
DR.  S.  THOMPSON 


2  ROCKWELL  INTERNATIONAL  SCIENCE 
CENTER 

ATTN:  DR.  S.  CHAKRAVARTHY 
DR.  S.  PALANISWAMY 
1049  CAMLSrO  DOS  RIOS 
THOUSAND  OAKS  CA  91360 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  luidertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  wiU  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-663 _ Date  of  Report  January  1995 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ _ _ 


4.  Specifically,  how  is  die  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ _ _ _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correcticm,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 

Organization 

OLD  Name 

ADDRESS  _ _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 
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